Quantum Simulation of Noncausal Kinematic Transformations 
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We propose the implementation of Galileo group symmetry operations or, in general, linear co- 
ordinate transformations, in a quantum simulator. With an appropriate encoding, unitary gates 
applied to our quantum system give rise to Galilean boosts or spatial and time parity operations 
in the simulated dynamics. This framework provides us with a flexible toolbox that enhances the 
versatility of quantum simulation theory, allowing the direct access to dynamical quantities that 
would otherwise require full tomography. Furthermore, this method enables the study of noncausal 
kinematics and phenomena beyond special relativity in a quantum controllable system. 
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Quantum simulations consist in the intentional repro- 
duction of a quantum dynamics on another quantum 
system that is, typically, more controllable [TJ. They 
promise to revolutionize computing technologies allow- 
ing us to solve otherwise intractable problems with min- 
imal experimental resources [2]. Several physical mod- 
els have already been proposed for quantum simula- 
tions: quantum phase transitions [3], spin models |U- 
[9] , quantum chemistry [10l [IT] , particle statistics includ- 
ing anyons [HI 03], many-body systems with Rydberg 
atoms [JJ], quantum relativistic systems [T5H25] . inter- 
acting fermion and fermion-boson [331 HE] models, 
Majorana fermions [271 [25] . the quantum Rabi model 
in superconducting qubits [29j . and relativistic quantum 
mechanics in circuit QED [30j . It is known that the com- 
putational power of a quantum simulator may overcome 
that of classical computers. However, the set of opera- 
tions that we can apply in the former is restricted com- 
pared with the versatility of the latter. For example, a 
wide set of unphysical but computable operations are not 
accessible to conventional quantum simulators. 

A quantum simulation can be seen as a process in 
which a quantum system is forced to behave according 
to a given mathematical model, closely reproducing its 
dynamics. At the same time, the simulator has a dynam- 
ics governed by the fundamental laws of nature. In this 
sense, we wonder whether quantum simulators may en- 
code processes violating their internal operating rules. In 
Ref. [31] . we gave a first example showing how to imple- 
ment quantum simulations of phenomena beyond quan- 
tum physics. Consequently, a natural question arises: is 
it possible to simulate processes violating special relativ- 
ity in a quantum device respecting it? 

In this Letter, we propose a formalism that allows the 
implementation of Galilean boosts and, in general, coor- 
dinate transformations, as spatial or time parity opera- 
tions, at any evolution time of a quantum simulation. 
This is significant to increase the versatility of quan- 
tum simulators, enabling the change of reference frame 
in situ during an experiment. The ability of generating 



these computable operations allows us to obtain corre- 
lations between different reference frames. These corre- 
lations include, among others, relevant physical quanti- 
ties as propagators and self-correlation functions, that 
would otherwise require full state tomography. More- 
over, one could also test and analyze the ultimate limits 
of a quantum simulator, exploring the exciting possibility 
of implementing noncausal kinematics as, e.g., instanta- 
neous translations or boosts in a controllable quantum 
platform. We will present the proposed method in the 
context of linear coordinate transformations, where the 
Galileo group is included. Moreover, this proposal also 
allows the implementation of nonlinear coordinate trans- 
formations, as is the case of accelerated frames. 

Let us first consider a linear transformation be- 
tween spacctime coordinates (t,x) — > (t',x'), i.e., x[ — 
Ylj=o a ij x i w ith (xo,xi) = (t,x), and the condition 
a oi = 0, assuring that t' is not mixed with x. Here, 
the coefficients a^ determine the kind of transformation 
that the quantum simulator will be able to reproduce at 
any evolution time of the dynamics. For instance, this 
could be a spatial (t, x) — > {t, ax) or time (t, x) — » {at, x) 
dilation, coinciding with spatial or time parity operations 
for a — —1. Another example may be a Galilean boost 
(t, x) — > (t, x — vt), where v is the relative velocity be- 
tween two different reference frames. We point out that 
the broad set of transformations with aoi = repre- 
sents the most general linear coordinate mappings im- 
plementable in a quantum simulator without prior knowl- 
edge of the evolved dynamics, as we explain below. 

To illustrate the formalism, we consider now a basic dy- 
namics governed by the equation idtip = —id x ij), where 
we assume natural units (h = c = I). This is a 1 + 1 
Dirac equation for a massless particle where, for simplic- 
ity, we have traced out the internal degrees of freedom. 
This dynamics takes place in the Hilbert space C 2 that 
we will call the simulated space. Our aim is to realize the 
above mentioned coordinate transformation during the 
unitary evolution, in such a way that our quantum state 
?p{x,t) is mapped instantaneously onto ^{x' ,t'). Notice 
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that this transformation is noncausal and, in general, vi- 
olates special relativity. Nevertheless, we will show that 
we are able to implement this unphysical mapping in an 
enhanced quantum simulator. 

The possibility of implementing this kind of operations 
in a quantum simulator arises from the fact that any wave 
function ip and operator 9 can always be expressed as 



and also 



ip(x,t) 



1 



(x,t)+^{x',t / )] + [iP(x,t) - i){x' ,t')] 



(x, t) = ~\ [9(x, t) + 9(x', t')] + [9(x, t) - 9(x', t')] 



Correspondingly, for the particular case 9 — d tiX , the 
time and spatial derivative operators are dt, x — \[dt,x + 
9f,x'] + \[dt,x — dt',x']- With these mappings, we can 
write the dynamical equation, idtip — — id x ?p, in terms 
of its even (e) and odd (o) components as follows, 

+ d° t ){r + r) = -m + d°)(r + r), m 

where ^° = \ ty(x, t)±i>(x', f)] , d\% = \[d t>x ± d t ,, x ,]. 
Now, we define a spinor ^(x, t) which belongs to C 2 £g)£ 2 , 
that we call the enlarged space, according to ^ = 
(ip e , if)°) T , where T is the transpose operation. The 
spinor VP is related to ip through the expression ip(x, t) — 
(1,1)^. An action as ip{x,t) — > ip(x',t') is, in gen- 
eral, forbidden in the simulated space. For instance, 
the transformation (i, x) —¥ (t, —x) is an operation that 
produces an instantaneous translation and inversion of 
a wavepacket centered in an average value xq to an av- 
erage value — xq. In this sense, it would violate causal- 
ity and the no-signalling condition for large enough xq- 
For example, this action applied on the wave function 
tp(x, t) — exp[— (x— x ) 2 /^/ x ], where ^ x corresponds to the 
Gaussian width, produces exp[— {x + xo) 2 /"f x }- This in- 
stantaneous operation is forbidden in the simulated space 
by special relativity. Nevertheless, the spinor W in the 
enlarged space, corresponding to ^(x^f), is just a z ^, 
i.e., ip(x',t') — (1, l)cr 2 ^/(x, t). This means that a physi- 
cal action like a z , acting on the enlarged space, gives rise 
to a physically-forbidden action on the wave function in 
the simulated space. 

The dynamical equation for ^(x,t) can be obtained 
from Eq. ^ separating its even and odd components, 
giving rise to 



d? df 



Qe go 
^ x ^ x 

d° d% 



(3) 



We can write d?'° in terms of d t and d x as follows, 



[d t ±d v ] = -[d t ±d aoo t] 



± andt T cn w d x ) 



(4) 



d x° = \ [ d x ± dx'] = \ [d x ± d {ailX+aiot) ] 



d x 



1 



anaoo 



-{±a 00 d x ) 



(5) 



We can substitute these expressions in Eq. ^ in order 
to obtain a Schrodinger equation for 



idt® = + a 2 <T x ]d x $>, 



(6) 
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ailZ 2(a°^ ai ° ■ In order t0 compute the dy- 
namics associated with Eq. one just has to define 
the initial condition for \&, i.e., &(x, 0) = ^[tp(x,0) + 
il>(x'(x,0),t'(x,0)),i(>(x,0) - ^(x , (x,0),t , {x,0))] T . No- 
tice that this state is completely determined by the ini- 
tial condition tp(x, 0) of the dynamics in the simulated 
space only if t'(x,0) = 0. Otherwise, prior knowledge 
of the full spatial and time dependence of the wavefunc- 
tion ip(x,t) is needed to determine ^(x,0). This is a 
circular argument given that the time dependence of if) 
is what we want to compute in the quantum simula- 
tion. We point out that the kind of initial condition 
with t'(x,0) = does not include Lorentz boosts in the 
genuinely relativistic regime. This is due to the fact that, 
in the case of relativistic Lorentz boosts, t'(x,t) is a lin- 
ear superposition of x and t. This means that, in gen- 
eral, the condition with vanishing t' is not fulfilled. Ac- 
cordingly, the set of allowed linear transformations with- 
out prior knowledge of the full wavefunction is given by 
t 1 = ctoat,x' = aiot + aux, including among others the 
Galileo group. Moreover, extensions to nonlinear coor- 
dinate transformations, as t' = aoot k ,x' = f(x,t), with 
f(x, t) a certain function of x and t, can also be con- 
sidered. This is because they have valid initial condi- 
tions as well, ip(x'(x, 0),t?(0)) — ip(f(x,0),0), i.e., con- 
ditions that are univocally determined by the knowledge 
ofip(x,0). 

Equation Q includes the dynamics of ij){x,t) and 
i/j(x',t'), allowing us to obtain any expectation value of 
both dynamics relating it to the measurement of observ- 
ables in the enlarged space. This can be expressed as 

(0)*M = WOW = <*l(})o(l >i)l*) 

= {*\(I + a x )®0\y), (7) 

(0)^, t >) = {^PW) = (*\a z ( J ) O ( 1 . 1 ) "*l*> 
= <¥|(I-(7 X )®0|¥>, (8) 

where we use (x\ip) — ij){x,t), (x'\ip') — ip(x',t'), and 
(x\^) = V(x,t). 

With the same tools we can obtain information about 
correlations between ip(x,t) and ip(x' , t') in terms of mea- 
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surements of observables in the enlarged space, 



Evolution 



Evolution 



0(1 A)<tJV) 



(oww,*) = wow) = m » , 

= <*|( £ r,-»a- 1 ,)®0|*). (9) 

The application of our protocol to systems of many 
particles [35] will give rise to measurements of correla- 
tions otherwise impossible to perform. This is because 
the process for obtaining the correlations involves the re- 
alization of two different quantum simulations, one for ip 
and another one for ip' . Subsequently, we have to imple- 
ment a full quantum tomography of both cases and finally 
store all the data in a classical computer and calculate the 
correlations. In general, if we are dealing with dynamics 
involving many particles, realizing the full tomography is 
impossible, while extracting the useful information can 
be accomplished applying our proposed ideas. 

Our method also allows us to include nonlocal and non- 
causal operations directly into the dynamics. For exam- 
ple, the equation 

id t ip = [icr x p x Il x + ma z + a x V{x)] ip, (10) 

where H x ip(x,t) = ip(—x,t), contains the spatial par- 
ity operator T\ x in the kinetic term. This means that 
the operation (x,t) (—x,t), (ccoo, ctoi, &io, = 
(1, 0, 0, —1) is included in the dynamics of the simulator. 



Equation ( 10 ) cannot be implemented directly in an 



experimental setup. Nevertheless, through the mapping 
ip(x, t) * = (ip(x, t) + ip(-x, t),ip(x, t) - ip(-x, t)) T , 
encoding properly the action of Tl x on Pauli operators, 
we obtain its image in an enlarged space, 



p<jy (j x ~t~ ml <£> a z + V e I ® <t x + V°(j x 



Here, V(x) = V e (x) + V°(x) is decomposed in its cor- 
responding even (e) and odd (o) parts. We consider 
now the nonrelativistic regime m ^> (|(p)|, |V(x)|), 
and the case in which the potential is explicitly odd, 
V(x) = V°(x). Then, Eq. (10 1 can be written in an 



interaction picture with respect to the mass term, and 
for a parity eigenstate associated to the ± eigenvalue, we 
have 



id t ip 



2m 



+ <y°y ± (d x v°) 



TP- 



(11) 



This equation describes the causal behavior of a nonrela- 
tivistic particle under the influence of an even potential. 
This shows an intriguing relation between the dynamics 



of Eq. (10 1, restricted to the case of odd potentials, and 



the dynamics of its nonrelativistic limit given by Eq. (Ill 



The first one shows a noncausal behavior evolving under 
an odd potential, while the second one is a causal equa- 
tion where the effective potential that emerges is even. 



o 
o 



c2 

V / 



u. 



Transformatic 



FIG. 1. Scheme of gates acting on the active (coloured 
circles) and control (white circles) qubits, implementable in 
different quantum platforms. In this example, we depict a 
quantum simulation process in which two kinematic trans- 
formations have been performed in the middle of the simu- 
lated dynamics. The first and fourth boxes (starting from the 
right), are used to generate the dynamics given by an evo- 
lution operator Ut- The intermediate boxes represent local 
gates U c i,c2 acting on the control qubits. 



Another example of the transformartions we can 
include is the time parity, (t, x) — > (— t, x), 
(aoo, ocqi, a\o, an) = {— 1,0,0, 1). Time is a global pa- 
rameter common to all quantum particles, which is in- 
troduced in quantum mechanics through the Schrodinger 
equation id t il> = Hip. We include this operation in our 
formalism through the equality 

ip(x, t) = i j [ip(x, t) +i>(x, -t)] + [ip(x, t) - i(>(x, -t)] 

(12) 

Notice that ip(x, —t) can be expressed as H t ip(x, t) — 
\[i>{x, t)+yj(x, —t)} — ^[ip(x, t)—tp(x, —t)], where H t is the 
parity operator. Accordingly, performing this mapping 
amounts to introducing a minus sign in the second term 



of Eq. (12 1. The operation n t cannot be realized in the 



real world because it implies instantaneous travel from 
the future to the past for a set of particles. Nevertheless, 
it can be implemented using an enlarged space as follows. 
The corresponding Schrodinger equation for H ^ H(t) 
can be written as 



id t (ip e + ^°) =if e (V> e + V>°), 



(13) 



where yb e = \[^[x,t) + i/>(x, -i)}, tp° = l[^ x ,t) - 
i/j(x, —t)] and H e = H, given that the Hamiltonian does 
not depend on time. This equation can be separated into 
even and odd parts and expressed in an enlarged space, 



id. 



(14) 



Notice that dttp e is odd with respect to the transforma- 
tion t —> —t, and d t 4> is even. Here, the relation between 
spinors in the simulated and enlarged spaces is again 
if) = (1,1)* with <F = (ip e ,ip°) T . Now, at any time of 
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the evolution, we can apply a time parity in the enlarged 
space multiplying the state at time t by the gate a z . This 
is due to the fact that ip(x, —t) = (1, l)a z ^>(x, t). In this 
case, to prepare the initial state in the enlarged space is 
especially easy given that 



V(x,t = 0) 



ip(x,t= 0) 




®ip(x,t = 0). 



(15) 



This is because at t 



0, r = ![#r,0)+#c,0)] = 
ip(x, 0) and V° = \ [ip{x, 0) - ${x, 0)] = 0. 

The inclusion of the time parity operation allows us to 
encode the expectation value of the propagator e" lHt in 
two observables in the enlarged space. To accomplish this 
task, we first evolve the initial wavefunction ^(x,t = 0) 
under the Hamiltonian a x ®H e of Eq. ( 14 ), generating the 
state V(x,t) = ~(ip(x,t)+ip(x, -t),ip(x,i)-ip(x, -t)) T . 
Later, we apply e - lI ® H " A to ty(x,t), producing 



$(M) 



1 ( i/>(x, 2t) + rl>(x, 0) 
xp(x,2t) - tp(x,0) 



Thus, according to the following equivalence, 



(16) 



(e- i2tH ) = (^(x,0)\i>{x,2t)) = (*K ^j(l,l)|*> 

= <*i(*»+w)i*>> (17) 

the measurement of a z and a y in the enlarged space will 
provide us with the expectation value of the propagator. 
In some cases, it is possible to relate the expectation value 



of Eq. (17) with self-correlation functions. For example, 



for spin systems one can write 



-i2tH 



) = (v^,t = o)| e -^|v(s^ = o)> 



(18) 



(ip'{s h t 



-i2tH 



a J #'( Sl ,t = 0)>, 



where Si refers to the spin degrees of freedom, \ip'(si,t = 
0)) = aj\tp(si,t = 0)), and <jj corresponds to u x ,a y ,cr z 
for j — 1,2,3, or, in general, to any Hermitian linear 
combination of them. In cases in which {H, tTj} = 0, it 



is possible to write the last line of Eq. ( 18 ) as 



= Q)\a j {-t/2)a j {tl2)\^{a i ,t = 0)), (19) 

with <Tj{±t/2) = exp (=pf H)aj exp (±f H). This corre- 
sponds to a self-correlation function that would require 
in general full tomography to be computed. 

Our protocol can be generalized in order to include 
the possibility of performing several transformations in 
a quantum simulation on one or many particles. For 
instance, given a wave function ip(x,y), where x and y are 
independent coordinates that may represent the position 
of two particles in one dimension or the coordinates of 
one particle in two dimensions, we have ij)(x, y) = ipe,e + 



ipe,o + i>o,e + i>o,o, where 

= \{ v) + V(-*> y)] 

+(-iy [ip(x, -y) + (-l)ty(-s, -y)] |, (20) 

and i,j ={e = 0,o = 1}. We then consider the spinor 
* = (^e,e,^e,o,ipo,e^o,o) T , that is related to ip(x,y) 
through ijj(x,y) = (1,1,1,1)$. The decomposition of 
Eq. p0| allows to generate easily operations of spatial 
parity in the x and y coordinates, W. x ip{x, y) = ip(—x, y), 
Uyip(x,y) = ip(x,-y), n x Ilyip(x,y) = tp(-x,-y). This 
is achieved applying just local gates to the state $ in 
the enlarged space, H x ip(x, y) = (1, 1, 1, 1) (cr z ® /) $, 
Uyili(x,y) = (1,1,1,1) (I ®a s )% U x U y i){x,y) = 
(1, 1, 1, 1) (<j z (8 a z ) $. In general, each inclusion of a 
new symmetry transformation in the quantum simula- 
tion amounts to doubling the Hilbert space. In different 
quantum optical implementations, like trapped-ion se- 
tups [53] and superconducting qubits [33], this implies 
adding another qubit encoding the proposed symmetry, 
see Fig. [TJ The interactions that appear due to the in- 
clusion of these symmetries are efficiently implementable 
in a digital quantum simulator with recently developed 
techniques [2~4"II55"] . 

In summary, we have explored the limits of quantum 
simulations via a formalism performing linear coordinate 
transformations during a simulated quantum dynamics. 
Among other features, we may compute spin temporal 
correlation functions without performing full tomogra- 
phy. Moreover, our method allows us to measure directly 
correlations between different reference frames. Finally, 
we may also study noncausal kinematics in a system re- 
specting the laws of quantum physics and relativity. 
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